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ABSTRACT 

We compute the complete gauge and chiral superheavy mass spectrum and cou- 
plings of the Minimal Susy GUT (based on the 210 — 126 — 126 — 10 irreps as the 
Higgs system) by decomposing SO(10) labels in terms of Pati Salam subgroup labels. 
The spectra are sensitive functions of the single complex parameter that controls 
MSGUT symmetry breaking. We scan for the dependence of the threshold correc- 
tions to the Weinberg angle and Unification scale as functions of this parameter. We 
find that for generic values of the GUT scale parameters the modifications are within 
10% of the one loop values and can be much smaller for significant regions of the 
parameter space. This shows that contrary to longstanding conjectures, high preci- 
sion calculations are not futile but rather necessary and feasible in the MSGUT. The 
couplings of the matter supermultiplets are made explicit and used to identify the 
channels for exotic {AB ^ 0) processes and to write down the associated bare d = 5 
operators (some of both are novel). The mass formulae for all matter fermions are 
derived. This sets the stage for a comprehensive RG based phenomenological analysis 
of the MSGUT. 



1 Introduction 

The Supersymmetric SO(10) GUT based on the 126, 126, 210 Higgs multiplets jUEl 
ElE] has, of late, enjoyed a much delayed bloom of interest motivated by its economy 
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and predictivity. Besides the traditional virtues of SO(10) this is the minimal renor- 
malizable model which has shown itself capable of matching the observed fermion 
spectra, including the prima facie GUT repellent feature of maximal mixing in the 
neutrino sector ISHOlCj- Beyond the traditional scenario of perturbative unification 
of couplings due to the RG flow between M s and M x it also offers strong indications 
that the gauge coupling becomes strong above the GUT scale. We have argued [HI E] 
that this necessarily leads to dynamical symmetry breaking of the GUT symmetry at 
a scale K v (just above the perturbative unification scale M v ~ 10 m GeV). Utilizing 
the quasi-exact supersymmetry at the GUT scale we made plausible [9 a scenario in 
which A[/ is calculably determined by only the low energy data and structural features 
of the theory (such as the gauge symmetry group, supersymmetry and the very re- 
stricted Higgs multiplets available to generate fermion masses -particularly neutrino 
masses- in a renormalizable theory). This scenario offers interesting possibilities of a 
novel picture of elementarity and dual unification characterized by a new fundamental 
length scale ~ A^ 1 characterizing the "hearts of quarks." [HI E]. 

The MSGUT is thus the focus of multi-faceted interest and a detailed phenomeno- 
logical analysis of the theory in terms of the structure dictated by its GUT scale 
minimality is thus called for. However such an analysis has been delayed by the 
computational difficulty of obtaining the GUT scale spectra and couplings and the 
effective Lagrangian describing the normal and exotic features (baryon and lepton 
violation etc) of the GUT derived MSSM (i.e extended by the leading (d = 5) ex- 
otic operators of the theory). The spectra and couplings are necessary to analyse 
threshold corrections to the gauge couplings near the GUT scale and are also a cru- 
cial input into deriving the lagrangian for exotic processes and parameters mediated 
by GUT scale massive fermions. In ^l] we presented techniques for computing the 
decomposition of SO(10) invariants in terms of the unitary labels of its maximal (Pati- 
Salam) sub group SU (4) x SU {2) L x SU (2) R . Once this decomposition is performed 
the computation of the complete spectrum and couplings is quite easy and the long 
standing vagueness regarding the "Clebsches" that arise can finally be banished. This 
allowed us to present, by way of illustration of the power of our method, the two most 
important mass matrices (4x4 and 5x5 respectively) affecting Electro weak sym- 
metry breaking, fermion masses and nucleon decay : namely those for the MSSM 
type Higgs SU(2)l doublets and baryon number violation mediating SU(3) C colour 
triplets that mix with the the doublets and triplets of the fermion mass (FM) Higgs 
(10, 126). Moreover since our methods allow computation of the actual couplings of 
Higgs to spinors we could also obtain the d = 5 operators for Baryon violation gen- 
erated via exchange of triplet Higgsinos contained not only in the traditional (6, 1, 1) 
submultiplets (of the 10 or those in the 126 [T2j) which had been noticed to provide 
a connection between neutrino masses and proton decay, but also in other channels 
arising from the exchange of colour triplets contained in (10, 1, 3)^26 submultiplets in- 
volved in neutrino mass generational]. A more complete calculation of these spectra 
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and effective lagrangians and an initial estimate of their effects is the subject of this 
paper. 

While the calculations presented were in their final stages we were collaborating 
and cross checking with another parallel calculation jTSj of mass spectra using a differ- 
ent [10J method which has since been published . Moreover another group [T^ ITo] has 
also recently published a calculation (using the same methods as d3J of spectra and 
baryon decay effective potentials recently. As far as computation of chiral superfield 
spectra are concerned our results coincide (upto normalization and phase conventions) 
with those of [T^]. However both our results diverge [TH ITB"] in certain details from the 
chiral spectra given in [T3j. Moreover as already noted by us in (an update to) [TT] 
we also disagreed with the results of ^3] regarding the Higgsino channels available for 
baryon decay in this model. We found [TT] that [JJ] obtained couplings between the 
126 multiplet and matter in the spinorial 16 representation which were in contradic- 
tion with ours^TJ not only as regards the numerical coeffcients but also in the heavy 
Higgs channels to which matter fields couple in a baryon number violating way. In 
the revised version of ^3] i.e this defect has apparently been corrected at least 
modulo disagreement on values of clebsches. We shall try to settle these questions 
by tracing the reasons for the continuing discrepancy in explicit detail and confirm 
our previous assertions. We have also analyzed the gauge Dirac multiplet structure 
arising from the super-Higgs effect and the masses and vevs responsible for the Type 
I and Type II mechanisms [TBI ITT] of neutrino mass generation. 

We emphasize that our method allows computation, not only of spectra but also 
of the couplings of all the multiplets in the theory (whether they are renormalizable 
or heavy-exchange induced effective couplings) without any ambiguity. Moreover 
our results are obtained by an analytic tensorial reprocessing of labels of fields in 
the Lagrangian. This approach might thus find preferment with field theorists in 
comparison with the more restricted capabilities of the approach of ^U], which, so 
far, has not proved capable of generating all the Clebsches of the SO (10) theory 
and which relies on an explict multiplet representative and computer based approach 
which is tedious to connect to the unitary group tensor methods so familiar to particle 
theorists. 

It has long been held by some that 5*0(10) GUTs specially Susy 5*0(10) GUTs, 
are essentially self-contradictory ^H] due to the apparently enormous threshold ef- 
fects that might arise due to the large number of superheavy Higgs residuals in these 
theories. Thus the authors of ^H] speak of the the "futility of high-precision calcula- 
tions in SO(10)". However these assertions have never been tested against any actual 
computations of mass spectra of a Susy SO(10) GUT and are only worst case esti- 
mates assuming that no cancellations occur. However we expect that cancellations 
will generically occur since the lepto-quark mass has no reason to lie at the edge of the 
mass spectrum nor are the coefficients all of the same sign. With the computed spec- 
tra now available, the threshold effects on observable quantities such as Sin 2 9 w , Mx 
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etc become computable in terms of the relatively small number [3] of GUT scale pa- 
rameters of the MSGUT. In fact [3] a single complex parameter controls the solutions 
of the cubic equation in terms of which all the superheavy vevs are defined. We have 
performed a preliminary scan of the parameter space of the MSGUT to see what is 
the typical size of these corrections. We find the striking result that such threshold 
corrections are generically < 10% of the 1-loop results [2H1 123 E3 EH] that underpin 
the the GUT scenario's viability. Moreover, for significant and possibly interestingly 
restricted regions of parameter space these corrections can be much smaller i.e as 
small as .5% of the one loop results. Thus far from indicating futility our results 
indicate that a thoroughgoing investigation of the compatibility of low energy pre- 
cision data with the threshold corrections may significantly constrain the parameter 
space of the MSGUT. In any case we show that inclusion of threshold corrections is 
necessary and not futile. Such an analysis is now in progress [27] . 

In Section 2. we present a brief review of the principal features of the minimal 
Susy SO(10) theory [TJ El El El and compute the gauge supermultiplet masses. In Sec- 
tion 3. we provide the PS reduction of the SO(10) Higgs superpotential. From this we 
computed the Chiral fermion mass terms and thus the supermultiplet spectra which 
we discuss here and list in an Appendix . In Section 4. we compute the threshold cor- 
rections to the 1-loop values of the unification scale Mx and sin 2 6w{Ms). In Section 
5. we present the couplings of the matter fields to FM Higgs fields in the superpo- 
tential as well as their couplings to the gauginos of the SO(10) model. This permits 
us to identify the possible channels for baryon violation in the low energy theory via 
exchange of Higgsinos or gauginos and compute the relevant effective lagrangians. 
Using the associated mass matrices we write down the d = 5 effective lagrangians for 
baryon and lepton number violation which arise via exchange of superheavy fermions. 
In Section 6. we discuss the mass formulae for the matter fermions in this model. The 
majorana mass terms of the left and right handed neutrinos and the £77(2) l triplet 
micro -vev responsible the Type II mechanism for neutrino mass is calculated along 
with the charged fermion mass matrices. In a final section we discuss our conclusions 
and results and plans for further investigations using the results derived here. 

2 The Minimal Susy GUT 

In accordance with our basic rationale we shall deal with a renormalizable globally su- 
persymmetric £0(10) GUT whose chiral supermultiplets consist of "adjoint multiplet 
type" (or AM) totally antisymmetric tensors : 210(<& iiJfci ), 126(£ iifeim ), 126(S ijfcZm ) 
= 1...10) which serve to break the GUT symmetry to the SM, together with 
Fermion mass (FM) Higgs lO-plet(Hj). The 126 plays a dual or AM-FM role since 
besides enabling Susy preserving GUT symmetry breaking, it also enables the gen- 
eration of realistic charged fermion masses and neutrino masses and mixings (via 
the Type I and/or Type II mechanisms) ; three spinorial 16-plets ^a(^4 = 1,2,3) 
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contain the matter supermultiplets together with the three conjugate neutrinos (vf;)- 
The 126 (£), 126(E) pair is required to be present together to preserve Susy while 
breaking U{1)r x U(1)b-l ~^ U{l)y and is capable [5J El U\ of generating realistic 
neutrino masses and mixings via the type I or type II seesaw mechanisms |TK| IT7j . 
The complete superpotential in this theory is the sum of 



iW — —MnHf + -TT&ijkl&ijkl + -^&ijkl&klmn&mnij + -jrT-^ijklm^ijki a. 



and 



+ ~T7 < ^ijkl^'ijmno^'klm,no + 'TrHi^jklmi'J^'ijklm + ^^ijklm) (1) 



W FM = h! AB i, T A cf 1 ^ B R l + ^f AB ^Ci% 1 ... li ^ B E il ... i5 (2) 



Our notations and conventions for spinors can be found in [TT] . The Yukawa couplings 
h' AB , f' AB are complex symmetric 3x3 matrices and one of them, say /' can thus be 
diagonalized (by an orthogonal transform UfU T using a Unitary matrix U which 
leaves the matter kinetic terms invariant) to a real positive diagonal form f' AB = 
F' a Sab, thus leaving 15 residual real parameters in Wfm- In addition the 7 complex 
parameters in Wam can be reduced to 10 real ones by absorbing 4 phases by Higgs field 
redefinitions. Then together with the gauge coupling one has in all exactly 26 non-soft 
parameters p]. Coincidentally, MSSM also has 26 non-soft couplings consisting of 9 
quark and charged lepton masses, 3 majorana neutrino masses 3 quark (CKM) and 
3 lepton (PMNS) mixing angles and 1 quark but 3 lepton CP phases together with 
3 gauge couplings and a /i parameter. Thus we see that the 15 parameters of Wfm 
must be essentially responsible for the 22 parameters describing fermion masses and 
mixings in the MSSM. 

The kinetic terms are given by covariantizing in the standard way the global 
SO (10) invariant D-terms 

1 — * — 1 

[■^-^{^ijklm^ijklrn + ^ij Um T>ijklm) + -^^ijkl^ijkl + H* Hi + ^* A ^ A ] D (3) 

Note that the extra factor of (1/2) achieves canonical normalization for the 126 
independent component fields of the self dual ( anti self-dual ) 126(126) representa- 
tions which would be otherwise be overcounted. We thus, unfortunately, differ from 
the normalization used in the parallel computations of JH] with which our results are 
nevertheless in agreement after appropriate rescalings and rephasings of parameters 
and fields. We emphasize that all our redefinitions of labels are unitary and thus 
maintain unit norm relative to the above kinetic terms. 1 



1 The relations between the quantities of (denoted by primes) and ours are S' = S/\/2, £ = 
S/\/2 (also for vevs) 7' = V2j, 7' = V27, rj = 2rj, M' = 2M. 
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The economy of the above superpotential eqns. (jll2|) is remarkable 4j. It's few 
couplings together with the functional flexibility of the chosen Higgs multiplet set 
and its ability (in common with other renormalizable models using just 10 — 126 FM 
Higgs) to fit the all the fermion mass datajHllZj, justify its claim to being the MSGUT. 
The "small" number of non-soft parameters (26 as in the MSSM) implies that after 
fitting Ej the known quark, charged lepton masses and quark mixing angles 
together with the neutrino mass splittings very little play is left in the model and 
it becomes predictive and thus falsifiable. The nearest related model(NMSGUT?)(in 
some ways more logically complete since all the FM channels allowed by renormaliz- 
ability would then be utilized) might be considered to be the one obtained by adding 
a 120-plet SO(10) FM Higgs. Alternatively one may consider SU(5) supplemented 
with right handed neutrinos or non-renormalizable terms j3]. Both models are are 
far less economical and not so predictive. Therefore, as advocated in detail in 
the first priority should be to pin down the predictions of this model. We began the 
development of a detailed framework for handling the group theoretic complexity of 
susy SO (10) models generally in JT] and this paper presents the results of calcula- 
tions using the techniques developed there for computing couplings and spectra for 
MSSM fields from the MSGUT tree action by decomposing the fields according to 
the SU(4) x SU(2) L x SU(2) R or Pati-Salam (PS) maximal subgroup. 

We now specify how the symmetry is broken down to the MSSM gauge group 
PU 121 E] by superlarge vevs contained in the 210, 126, 126-plet scalar vevs. Before 
doing so we introduce our submultiplet naming and indexing conventions. A host of 
further details related to the Pati Salam decomposition of SO(10) can be found in 
our earlier paper JT] where the foundation for the current program of computation 
of states, masses and couplings of this theory was laid and the spectrum of MSSM 
like SU(2)l doublets and SU(3) C "baryon decay" triplets first computed. 

We denote quantum numbers w.r.t the SM gauge group by enclosing them in 
square brackets while those with respect to the PS group are denoted by round 
brackets. We have adopted the rule that any PS submultiplet of an SO(10) field 
is always denoted by the same symbol as its parent field, its identity being estab- 
lished by the indices it carries or by additional sub/superscripts ((a), (s), ±, L, R) 
denoting (anti-)symmetry or (anti)-self duality, if necessary. On the other hand, 
since one often encounters several chiral MSSM multiplets of the same type aris- 
ing from different SO(10) Higgs multiplets we will also introduce a naming con- 
vention using roman letters for these multiplets. If we need to denote the scalar 
component of a chiral superfield we use a tilde over the superfield symbol and some- 
times use a superscript "F" to denote fermionic components of chiral superfields 
while gauginos are denoted by A. Our notation for indices is as follows : The real 
indices of the vector representation of SO(10) are denoted by i,j = 1..10. The 
real vector index of the upper left block embedding (i.e. the embedding specified 
by the breakup of the vector multiplet 10 = 6 + 4) of SO(6) in SO(10) are de- 
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noted a,b — 1,2. .6 and of the lower right block embedding of SO (4) in SO(10) by 
a, $ — 7, 8, 9, 10. These indices are complexified via a Unitary transformation and de- 
noted by a,b = 1,2,3,4,5,6 = ~p, ~p* = 1, 1*, 2, 2*, 3, 3* where 1 = 1,2 = 1* etc. Simi- 
larly we denote the complexified versions of a, (3 by a, (3 = 7, 8, 9, 10 = 0. Using this 
complexification we showedjlT] how all SO(6) x SO(4) subinvariants of SO(10) tensor 
invariants could be systematically converted to SU (4) x SU (2) L x SU(2) R invariants 
whose indices are as follows : The indices of the doublet of SU(2)l(SU(2) r ) are de- 
noted a, j3 = 1, 2(d, P = 1,2). Finally the index of the fundamental 4-plet of SU(4) is 
denoted by a (lower) fM, u — 1, 2, 3, 4 and its upper-left block SU(3) subgroup indices 
are /i, v = 1, 2, 3. The corresponding indices on the 4 are carried as superscripts. 
These doublets and quartets correspond to the chiral spinor representations of the 
SO (4) and SO (6) subgroups of SO (10). Details of the spinorial invariant decom- 
position techniques may be found in Ilj. The component of the SU(4) adjoint in 
the direction of the Gell-Mann generator is labelled with a superscript (15) or 
(B-L). 

Thus the PS decomposition of our SO(10) multiplets is 

= 210 = 0/(15, 1, 1) + 0(1, 1,1)+0^ L) (15, 3,1) +0^(15, 1,3) 

+ <Wd(a)(6,2,2) +0^ (S) (1O,2,2) + 0% S) (TO,2,2) (4) 

S = £+ = 126 = S^ (s) (T0, 1,3) + £^ (s) (10, 3, 1) + Sg(6, 1, 1) + £^(15, 2, 2) (5) 

£ = E" = 126 = E^! S) (10, 1, 3) + C (S) (T0, 3, 1) + EJJ (6, 1, 1) + ^(15, 2, 2) (6) 

H = 10 = tf ad (l, 2, 2) + tf M „(6, 1, 1) (7) 
q A = 16 = 16+ = (4 + , 2+) + (4_, 2_) = F£(4, 1, 2) + F^(A, 2, 1) (8) 



F(4, 2, 1) = (Q Aa , L a ) F(4, 1, 2) = (Ql, L a ) (9) 

with 




The GUT scale vevs that break the gauge symmetry down to the SM symmetry 
are 012]: 

((15, 1, 1)) 21Q : (<f)abcd) = ^abcdeftef (H) 

where [e e f] = Diag(e 2 , e 2 , e 2 ), e 2 = ir 2 . One dualizes 

0ab = T7 tabcdef 4> 'cdef (12) 
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Then in SU(4) notation [<p u ] this vev is 

[(tf)] = *Diag(I 3 -3) = ^ (13) 

((15, l,3)) 2 i : i&ab&p) = (14) 
where [e-s] = Diag(e 2 , e 2 ) which translates to 

((1, 1, l)) 2 io : ((f> & pjs} = P e &fa~s (I 6 ) 



((10, 1, 3)) m : (E 13g§6 ) = a = -i<E$ +) ) = ^ (17) 



S 44 

((TO, 1, 3)) 126 : (Z mn ) = a = ? (S^ 44 ) = - J 
Substituting these vevs into the superpotential one obtains 



W = m(p 2 + 3a 2 + 6u 2 ) + 2A(a 3 + 3pu 2 + 6au 2 ) 

+ (M + r/(p + 3a-6cj))(ra (19) 

the nontrivial F term conditions are thus : 

2mp + 6Xu 2 + r]aa = (20) 

2ma + 2\(a 2 + 2u 2 ) + r/aa = (21) 

2muj + 2uj\(p + 2a) - r]aa = (22) 

(M + r](p + 3a - 6u))a = (23) 

The vanishing of the D-terms of the SO (10) gauge sector potential imposes only 
the condition 

\a\ = \a\ (24) 

Except for degenerate cases corresponding to enhanced unbroken symmetry 
(SU(5) x U(l), SU(5), (j3,2,2,b-L) G3,2,r,b-l etc) [U [J3] this system of equations is 
essentially cubic and can be reduced to the single cubic equation jlj for a variable 
x = —\uj/m : 

8x 3 - 15x 2 + Ux - 3 = -f (1 - x) 2 (25) 
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where £ = — and the other vevs can be expressed in terms of values of the variable 
x which solve eqn ()25|) .This parametrization of the MSGUT ssb problem [3] is of great 
help computationally and clearly exhibits the crucial importance of the £ i.e of the 
ratio M/m. The important role played by a similar ratio in the other renormalizable 
SO(10) GUT based on the 45, 54, 126, 126 representations has already been noted 

When we measure vevs or masses in units of m = m/X we will put a tilde over 
the symbol. We also define the additional dimensionless parameters fj = rj/X and 
M H = M H /m. 

Then the dimensionless vevs are u = —x and 



{x 2 + 2x-l) „_x(5x 2 -l) 2x{l -3.x)(l + a; 2 ) 



[1 — x) (1 — X) 2 fj (1 — x] 



The solutions of the cubic equation (|25|) are generically complex. We will therefore 
nowhere assume hermiticity for our mass matrices, preferring to leave them undiago- 
nalized for eventual numerical diagonalization so that all our results are applicable in 
the general case. We will not generate arrays of expressions in terms of the variable 
x , although it is easy to do so since, practically speaking, the substitutions are now 
handled via a computer anyway. 

We conclude this section with a description of the super-Higgs effect for the break- 
ing SO (10) — > SU(3) x SU(2) x U(l)y which is achieved by the above superheavy 
vevs. As is well known, as a consequence of gauge symmetry breaking, each massive 
gauge boson forms a massive supermultiplet together with its longitudinal goldstone 
pseudo scalar (and its real scalar partner) as the 4 bosonic degrees of freedom. It's 
gaugino and the chiral fermion superpartner of the Goldstone scalar pair make up 
one Dirac fermion super-partner also with 4 degrees of freedom. This is the so called 
Dirac or massive vector gauge supermultiplet. These gauge boson/gaugino masses 
are the most fundamental thresholds of the GUT and it is appropriate to begin with 
a discussion of their values for this model. In Section 4. we follow |2~T] E2j to compute 
the threshold corrections using the spectra we compute. Then a Dirac gauge coset 
multiplet in representation R of the MSSM gives rise to a RG mass threshold above 
which the gauge and chiral components of the Dirac multiplets separately contribute 
—3Ti(R) and Ti(R), respectively, to the beta function coefficients of the individual 
MSSM couplings (including the U(1)y coupling !). 

It is easiest to keep track of the gaugino masses and mixings. The combination 
of chiral fermions that forms a Dirac fermion together with a gaugino must, for 
consistency, be a zero mode of the mass matrix arising from the superpotential and 
this makes it easy to disentangle the gauge spectrum even in the case of complex vevs 
and parameters. For the symmetry breaking to the MSSM the gauginos of the coset 
SO(10)/Gs2i lying in the PS representations (6,2,2) © (1, 1,3) plus the triplets and 
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anti-triplets in (15, 1, 1) (i.e 33 Dirac multiplets in all) obtain a mass by pairing with 
chiral AM Higgs fermions. One need only substitute the vevs given above into the 
PS decomposition of the gaugino Yukawa terms which have the form 

9^{^ < ^ijkl > ^jm^mjfc/ + 7^j'( < ^ijkln > ^im^mjkln+ < ^ijkln > ^im^mjkln) }+H.C 

(27) 

One finds the following gaugino masses : 

• (i) G[1,1,0] : m Xa = VlOg\a\. 
The mass term is 



3 (v /2A(«°) - V3\^)(a*^_+a*i: 4m+ ) + H.c 
v2 

= m XG 9±(e-^G 4 + e-^G 5 ) + H.c. 



Ge . (yf a( " o) - Vl A(15)) 



G. = ^ ; G .= ^- (28) 

The naming conventions for the chiral states are given in Section IV and the 
Appendix . Here 7 CT , 75- are the phases of a, a. Since the representation is real, 
the mass matrix Q in this sector is symmetric. The complete G[l, 1,0] sector 
mass matrix(includiing gauginos)^ is 6 x 6 while its pure chiral part G (which 
arises only from the superpotential) is 5 x 5 and symmetric and the 5-tuple 
(0, 0, 0, a, a) is both a left and right null eigenvector of G - as will be obvious 
when it is presented further on (Section 2. and Appendix). 



ii) J[3, 1, -§] © J[3, 1, §] : m Xj = g^8\a\ 2 + 1Q\uj\ 2 + 2\a\ 2 

In this case (J^p = A A (J4) ^ = A// 1 pair up with the combinations correspond- 
ing to the left and right null eigenvectors Vqjl = Nj(—a,2a,2\/2uj), v { 



OJR 



Nj(a, 2a, 2^2uj) of the complex, non symmetric, upper left 3x3 sub-matrix J 
of the 4x4 mass matrix J in the J sector. The gaugino mass terms are 

^J 4 (2v^aV 2 +4wV 3 -v^^Ji)-^(2y2a*J2+4^*J3 + v / 2(7*Ji)J 4 + J H".c (29) 

iii) F[l, 1,-2]© F[l, 1,2] ; m XF = g^2A\iu\ 2 + 2\a\ 2 

The chiral partners of the gauginos F 3 = \ R+ ,F 3 = A R _ correspond to the 
right and left null eigenvectors vqfr = (— cr, vT2iu;) T ; vqfl = (a,y/T2iu) of 
the 2 x 2 F — F chiral fermion mass matrix. The mass terms are 
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- gF 3 (iu* V2AF 2 + V2a* F 1 ) + g (ioo* \[2AF 2 - V2a* F ± ) F 3 (30) 

• iv) E[3, 2, -§] © E[3, 2, ±] ; m XE = gy J(4\a - cu\ 2 + 2\w - p\ 2 + 2\a\ 2 ). 

The chiral partners of the gauginos E 5 = X^ia, E 5 = correspond to the null 
eigenvectors v er = («cr, \/2(a — w), p) T ; v el = ( — 2<r, \/2(a — cj),c<j— p) 
of the upper left 3x3 corner E of the E sector 4x4 chiral fermion mass matrix 
S. -Ei, Ei do not mix with other E-sector multiplets. The mass terms are 

gE 5 (-iV2a*E 2 + 2(a* - u*)E 3 + V2(oo* - p*)E 4 ) 
+ g((iV2~W*E 2 + 2(a* -lu*)E 3 + V2(lu* -p*)E 4 )E 5 (31) 

• v) X\3, 2, §] © X[3, 2, -§] : m Xx = g^A\a + u\ 2 + 2\p + oo\ 2 . 

The chiral partners of the gauginos X 3 = Xn i2a ^ -^3 = correspond to the 
null eigenvectors Vqxr = ( _ v2(a + u>), uj + p) T = vqxl of the upper left 2x2 
corner X of the 3x3 X-sector chiral fermion mass matrix X . The X-gaugino 
mass terms are : 

gX 3 {-2(a* + o/)X! + V2(p* + u*)X 2 ) + g(-2{a* + u*)X 1 + V2(p* + u*)X 2 )X 3 

(32) 

3 AM Chiral masses via PS 

Our approach to opening up the maze of MSSM interactions coded in the deceptive 
simplicity of the SO(10) form of the GUT action is to rewrite SO(10) invariants as 
combinations of PS invariants using the translation techniques developed by us jTTj. 
Although tedious, our approach allows one to keep track of all phases and normal- 
izations without any ambiguity. Once this is done making contact with the MSSM 
phenomenology becomes trivial since the embedding SU(3) x SU{2)l x U(1)y C 
577(4) x SU(2)l x S77(2)_r is trivial and transparent if one keeps in mind that 

Y = 2T m + B - L 
We obtain for the PS form of the different terms in Wam 

^Im = -™{<p;w+<p«ts)^ ( 33 ) 
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-^■^ijkim^ijkim — M{E^ EjJ + 2E^ '"""S^ ad + .E^j ) + s|J (L) . E (s)(L) ) } (34) 

^A0 3 = A[-|0;0,>/-2.(0;0^ (s) A Q t ) ) 

- ^{0/(0^.0^+0^.0^)} 

+ {(0 ^u(R)a$ ^A(s)a^(L)a/3)) 

+ ^ + ^^A^))} 

+ v / 2{0^ ) (0^>; i?)Q/3 + 

— » , — * — » , — * 

- <t>{{<t>v{R)4\{R) ~ <t> v (L)4\\L))} 

+ ^{^(^AiA?R)7a + 0ML?^ A (2)/30A^) 7 a}] ( 35 ) 



= 7 [^ (a) r A(a V + ^;sr^ 

+ (-H ,lUW (() u{R) .T lfl x( S )(R) + ^( a )0A ! (L)-S (L) ) 

~ 2^ ^fW^»« 

- ^„(a)^ A(a)Qa S/ aa + ^#^ (a) E Ma) ] (36) 



l 7 0SiJ = 7[-^^(a) S ' A(a VA M + 0;Sr^ aa 

- -^{^ (a) s^0, A(s)aci + H^vr*^} 
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+ i^a^^^ + l^^} (37) 



|^EE = 77[2^(E^E A ; d + E^a G , , 



+ ^^(S^aw^-E^^ + E^ (s)(L) .E (s)(L) ) 

JZ* -» -» ->\1V 

+ 0(E^ {s)(R) .E(^ ){i?) - S^ (s)(L) • E (s)(L) ) 

=f rad/, /3 v mA(s) , v (s) ,uA(s)/3x 
Z22 > ^A( S )a^d/3(i*) + v\aj3(L)Qa ) 
~ -=i - _=ynA(L) 

2(E^^ (R) .E^ (i?) + E Ma) 0^ (L) .E ) 
+ 2(r Ka) ^( L ).S,A(L) + E, Ka) ^ (R) .E^^) 

- 2 v^r%^ )a , + e^e^0- (l)q/3 ) 

^/9^ a/3yr(s)(«)7 v ^A(s) , ^ a/3 v (s)(L) 7 ^A(s) s 

- V A</V ( R ) % A/ 3 ^(fl) 7 d + (L) %A/3 ^(L) 7 cJ 

_i_ ,- A /o/' w^^V A _i_ v v P^vKs) \ 



The purely chiral superheavy supermultiplet masses can be determined from these 
expressions simply by substituting in the AM Higgs vevs and breaking up the contri- 
butions according to MSSM labels. 

It is again easiest to keep track of Chiral fermion masses since all others follow 
using supersymmtery and the organization provided by the gauge super Higgs effect. 

There are three types of mass terms involving fermions from chiral supermultiplets 
in such models : (A) Unmixed Chiral (B) Mixed pure chiral (C) Mixed chiral and 
gaugino . 

3.1 Unmixed Chiral 

A pair of Chiral fermions transforming as SU (3) x SU (2)^ x U (l)y conjugates pairs 
up to form a massive Dirac fermion . For example for the properly normalized fields 

A[l,l,-4] = ^p A[l,l,4] = -££ (39) 
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one obtains the mass term 



2(M + r](p + 3a + 6u))AA = m A AA 

The physical Dirac fermion mass is then \mA\ since the phase can be absorbed by a 
field redefinition. By supersymmetry this mass is shared by a pair of complex scalar 
fields with the same quantum numbers. If the representation is real rather than 
complex one obtains an extra factor of 2 in the masses . There are in fact 19 types 
of such multiplets and their (roman letter) lables are given along with their masses 
and SO(10) origins in Table I in the Appendix. The case of the sectors C[8,2, ±1] 
and D[3, 2, ±1] bears special mention. The mass terms for these multiplets arise only 
between pairs drawn one each from £(15, 2, 2), £(15, 2, 2) and there is no mixing 
between a C, C or D, D drawn from the same SO(10) multiplet simply because the 
superpotential does not contain any term containing £ 2 or £ . This was the reason 
for the discrepancy in this sector between the results of [HI H3j and [Hj : there simply 
is no such mixing. 



3.2 Mixed Pure Chiral 

In this case there are no contributions from the gaugino Yukawas or the D-terms to 
the supermultiplet masses, but there is a mixing among several multiplets of the same 
SM quantum numbers. There are only three such multiplet types : 

. a) [8,1,0]^,^) = (</>/, <^ 0) ) 
These mix with mass matrix 



/ (m-Xa) -V2Xu \ 

\ -V2Xu m + X{p-a) J v 1 

with both rows and columns labelled by (i?!,^)- The masses are the magni- 
tudes of the eigenvalues of the matrix 1Z . 



\K ± \ = 2|m[l + (| - 5) ± y^) + 2<D 2 ]| = m R± (41) 
While the corresponding eigenvectors can be found by diagonalizing the matrix 

The mass matrices of the electroweak doublets h[l,2, 1], h[l,2, —1] and colour 
triplets t[3, 1, — |],f[3, 1, |] which mix with the multiplets contained in the 10 
plet FM Higgs are the most crucial ones for determining the phenomenology 
of the effective MSSM that arises from this GUT. These matrices were first 
calculated in [H](v2) and later, stimulated by a contradiction with a recent 
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paper the d — 5 baryon violating operators induced by the exchange of 
heavy Higgsinos were computed and added to a revised version JT] (v4) by using 
the Clebsches for the 16 • 16 • 126 and 16 ■ 16 • 10 invariants calculated earlier by 
us. Thus one has : 

• b) [l,2,-l](Ai,/l 2 ,/l3,/l4)©[l,2,l](/il,/l2,/l3,/i4) 
-(fra V (15)a yKl5)« 4>2i\ ffi(T T . y(15) v (15) <^4i 

These multiplets label the 4 rows and columns of the 4x4 mass matrix 7i|llj 
which is given in the collection of mixing matrices in Appendix I . We note that 
we have redefined our mass parameters m, M by a factor of 2 relative to those 
we used in [TT] . To achieve the MSSM spectrum of one pair of light doublets, it 
is necessary to fine tune one of the parameters of the superpotential (e.g M#) 
so that DetTi = 0. By extracting the null eigenvectors of H^H and H.T0 one 
can compute the coefficients of the various bi-doublets in the light doublet pair, 
and, in particular, we can find those for the doublets coming from the 10, 126 
multiplets which couple to the matter sector (see Section 6.). In this way the 
SO(10) constraints on the fitting of the Yukawa coupling matrices h' AB , f' AB 
can be brought into focus and the invalid assumption that the squares of these 
coefficients Ej add up to 1 can be dispensed with. 

• (c) [3,1, §](*!, i 2 ,h,U,is)® [3, l,-§](ii,i 2 ,t3,f4,f 5 ) 

For generic values of the couplings all these particles are superheavy. These 
triplets and antitriplets participate in baryon violating process since the ex- 
change of (t\, £ 2 , ti) © (ti, F 2 ) Higgsinos generates d = 5 operators of type QQQL 
and luud. The strength of the operator is controlled by the inverse of the t — t 
mass matrix T which we computed in [TI] and is given in the Appendix . We 
shall examine how d = 5 baryon and lepton number violating operators are 
generated in Section 5.. 

3.3 Mixed Chiral-Gauge 

Finally we come to the mixing matrices for the chiral modes that mix with the 
gauge particles as well as among themselves. Apart from threshold effects, these are 
of some interest since one might ask whether the new types of coset gauginos present 
in in SO(10) but not in SU(5) namely SO(10)/SU(5) ~ E[3, 2, |] © £[3,2,-§] © 
F[l, 1, 2] © F[l, 1^-2] © G[l, 1, 0] © J[3, if] © J[3, 1, -f] (the SU(5)/G 32 i leptoquarks 
are X[3, 2, — |]©X[3, 2, |] ) might not mediate interesting exotic processes by inducing 
d = 5 operators via mixed gaugino-chiral exchange. We have examined this question 
in some detail in Section 5. 
These multiplet sets are : 
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a) [1, 1, 0] (G 1 , G 2 , G 3 , G 4 , G B , G fl ) = (0, <P^\ 0g, ^5±> , ^ A(a y Al5) ) 

which mix via a 6 x 6 mass matrix given in the Appendix. The complex 
conjugates of the 6th row and column form left and right null eigenvectors 
vogl,vogr of the upper left 5x5 block G of Q. The determinant of Q is 
generically non zero although the determinant of the submatrix G vanishes. It 
will clearly not affect the evolution of the MSSM gauge couplings at one loop 
due to the singlet quantum numbers . 

b) [3, 2, -\\{E 2 , E 3 , E 4 , E 5 ) © [3, 2, \\{E 2 , E 3 , E A , E 5 ) = 

( S 4ai> S \ fia^' A S) ® ( S /a2' ^jLaV ^pLi' 

The 4x4 mass matrix £ ((Ei, E\) = (S^j, ^ 4 ^) do not mix with the others) has 
the usual superhiggs structure : complex conjugates of the 4th row and column 
are left and right null eigenvectors of the upper left 3x3 submatrix E. £ has 
non zero determinant although the determinant of E vanishes. As for the case of 
C[8, 2, ±1] and D[3, 2, ±7/3] type multiplets one finds that the conjugate types 
of E type multiplets drawn from the same SO(10) representation cannot mix. 
Furthermore explicit computation using the decomposition of the superpotential 
given in Section 3. shows that £>i[3, 2, ~] = S^ a2 and £a[(3, 2, — ~] = S^j in fact 
decouple from the other members of the E sector so that the E sector mixing 
matrix is 4 x 4 (including gauginos) and 3x3 excluding gauginos. Note that our 
assertion is not that these couplings cancel but simply that they do not appear. 
To see why, for instance, there is no term mixing say E\ = £^ q2 with E 3 = 0^,1 
coming from the $(10, 2, 2) we observe that the terms mixing S(15, 2, 2) and 
<&(decuplet, 2, 2) via a righthanded vev could only come from the following two 
terms in eqn.(|38|) : 



We see that the pairs £(10, 1, 3) and 0(10, 2, 2) and S(T0, 1, 3) and 0(10, 2, 2) 
simply do not mix. Now it is obvious that a £ right handed vev will mix only E 2 
coming from £(15, 2, 2) with E 3 coming from $(10, 2, 2) but not E\ coming from 
£(15, 2, 2) with E 3 coming from $(10, 2, 2). Similar considerations account for 
the other decouplings between E\,E\ and the rest of the E sector. Ultimately 
this correlation is accounted for by the correlation between the duality prop- 
erties of SO (6) decuplets and £0(4) triplets within the SO(10) self-dual and 
anti-self-dual multiplets £, £. 

• c) [1,1, -2]^, F 2 ,F 3 )® [1,1, 2](F 1 ,F 2 ,F 3 ) 
= (Em(ro), 0(ij-), \b-)) © 0(/e+)' \r+)) 
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The mixing matrix T has the usual structure . The residual massive eigenstates 
after separating off the the two Dirac fermions of mass 



m Xp =g(24\uj\ 2 + 2\a\ 2 )i (42) 

is a Dirac fermion of mass 

m F = \l\J\ a \2 + 121^2 (43) 
and the form of its chiral parts is 



F = N F (iVV2uoF l + aF 2 )e i ^-^ ) 
F = N F (-iVUuF 1 + aF 2 ) 



Np 1 = ^I2\u\ 2 + \a\ 2 (44) 

• d) [3, 1, -|](Ji, J 2 , J 3 , J 4 ) © [3, 1, |](Ji, J 2 , J 3 , J 4 ) 

= > , K) © (S A 4(R+) , 0/, 0^0) > V) ' 

The 4x4 mass matrix J has the usual super-Higgs structure : complex conju- 
gates of the 4th row and column are left and right null eigenvectors of the upper 
left 3x3 submatrix J. J has non zero determinant although the determinant 
of J vanishes. 

• e) [3,2,|](X 1 ,X 2 ,X3)©[3,2,-|](X 1 ,X 2 ,X 3 ) 

mix via a 3 x 3 symmetric matrix X so the left and right null eigenvectors of the 
upper left 2x2 submatrix X, formed by the complex conjugates of the third 
row and column of X, are the same. Separating off the two Dirac [3, 2, ±|] 
gauge fermions of mass 

m Xx = g\jA\a + u\ 2 + 2\p + u\ 2 
one is left with a Dirac fermion of mass 



2j 


m\{2\ 


x\ 


2 + l 


l-x\ 


2 ) 




|] 


— X 





m x = 2{2\m + X(a + u)\ + \m + Xcu\) 
whose chiral parts are also neatly expressed in terms of x 

(X,X) = 1 =(e^' m -^\V2xX 1 ^l-x)X 2 ),(V2xX 1 + (l-x)X 2 )) 

^J2\x\ 2 + |1 -x\ 2 

(45) 



17 



This concludes our description of the superheavy mass spectrum of the minimal 
susy GUT. As mentioned earlier our results were calculated in collaboration with 
the authors of^3j and are in agreement with the chiral spectra calculated in ^3] : 
whose results also confirm our earlier results ^I] on the phenomenologically important 
matrices 7i, T. Moreover we have evaluated the mixing of the gauginos with the chiral 
fermions explicitly and calculated the gauge spectra and eigenstates besides furnishing 
all the couplings in the superpotential sector explictly. The gauge couplings and the 
gauge Yukawa couplings to matter will be be given in the Section V. The gaugino 
mixing with the chiral fields will be useful to us when we examine B + L violation 
mediated by gaugino exchange as well by Higgsino triplet exchange and when one 
wishes to examine the flow of gauge couplings past the gauge thresholds. 

4 RG Analysis 

The first phenomenological success of GUTs was the 1-loop calculation of the numer- 
ical value of Weinberg angle [2H|. This was followed by the prediction [21] and then 
the verification [2S] of an amazingly exact compatibility between UV gauge coupling 
convergence in the MSSM and the precision LEP data. The large mass at which the 
top quark was eventually discovered and the associated large value Sin 2 6w ~ -23 
verified the originally somewhat far fetched conjecture of [21] : a historical fact that 
is still not always appreciated. The proposal of Weinberg [21] for calculating threshold 
effects within an effective field theory picture using mass independent renormalization 
schemes such as the standard MS renormalization scheme was taken up and developed 
in detail in [22]- Thereafter, using these results, it was argued [TH] that high-precision 
calculations in SO(10), and particularly in supersymmetric SO(10) models which used 
large representations such as 210,54,126 etc, were futile. This was due to the huge 
corrections to the one loop predictions that they expected in view of the large number 
of superheavy fields and the expected span in their masses. It should be remarked 
however that without an explicit calculation cancellations that might naturally occur 
would be overlooked. Such calculations were never done. These negative expectations 
were a motivation for the development |29| of a whole genre of SO(10) models that 
eschewed large representations (and thus parameter counting minimality) in favour of 
models with a plethora of small representations and non-renormalizable interactions. 

The other approach [TJ [21 El E3 EDI H] approach has all along been to retain 
renormalizability of the fundamental theory. We regard retention of Higgs multiplets 
just adequate to account for the gauge and fermion spectrum via renormalizable 
couplings as a sine qua non for even being clear as to what is testable about a given 
model. The inverse approach where representations ("hypotheses") are multiplied 
without necessity seems regressive to us. Thus the proposal of the Susy SO(10) 
GUT based on the 210 - 126 - 126 - 10|H [2] Higgs system as being the Minimal 
Susy GUT [3J E] must live or die by the criterion : Are the one loop values of 



18 



Sin 2 9w and Mx generically stable against superheavy threshold calculations ?. By 
"generically" we mean : for a non-singular subset of the parameter space. So far 
this question could not be answered definitively since no complete mass spectrum 
was available in any Susy SO(10) model to settle the issue. Partly this was due to 
the lack of accessible techniques to calculate mass spectra and couplings in these 
models due to the difficulty in obtaining the relevant SO(10) "Clebsches". Over the 
last few years we have developed a complete technology for translating SO(10) 
tensor and spinor labels into those of the unitary labels of the Pati-Salam maximal 
subgroup SU{4) x SU(2) L x SU(2) R of SO(10). This allowed us to compute first 
the mass matrices of SM type doublets and proton decay mediating triplets and then 
the complete spectrum and couplings reported in this paper. The partial technology 
of (HUE! nas a l so been used to compute [131 E] this spectrum (but not the couplings). 
With the correct spectrum in hand we can apply the standard formulae of Hall [22] to 
compute the changes in the 1-loop GUT predictions as functions of the few MSGUT 
parameters (£ = MX/mrj, A, rj, 7, 7, m, g, Mu) which are relevant at the GUT scale. 
Thereafter we can scan the parameter space to see how the corrections vary with 
these parameters. 

A few remarks on the role of the parameters are in order. The parameter £ = 
XM/rjm is the only numerical parameter that enters into the cubic equation eqn. (|25j) 
that determines the parameter x in terms of which all the superheavy vevs are given. 
It is thus the most crucial determinant of the mass spectrum. The dependence of 
the threshold corrections on the parameters A, rj, 7, 7 seems quite mild (logarithmic) 
(this is especially obvious for the unmixed chiral multiplets) and is also suggested 
by our preliminary scans of the parameter space. Thus changing A, 77 by a factor 
of 100 each yields plots vs £ that seem indistinguishable from the ones presented 
below. From equations (j2(ij) we see that m/A can be extracted as the overall scale 
parameter of the vevs. Since the threshold corrections we calculate are dependent 
only on (logarithms of) ratios of masses the parameter m does not play any crucial 
role in our scan of the parameter space : it is simply fixed in terms of the (threshold 
and two loop corrected) mass My = Mx of the lightest superheavy vector particles 
mediating proton decay : which mass is chosen, in the approach of Hall, as the 
common "physical" matching point in the equations relating the running MSSM 
couplings to the SO(10) coupling[22j. Inasmuch as we take the parameters A, 77 as 
given, and the parameter m is set by the overall mass scale, the freedom in the 
parameter £ is essentially that of choosing the 126 — 126 mass parameter M i.e the 
freedom in choosing the dimensionless parameter £, is essentially that of the ratio 
M/m : which ratio is already known to be a crucial control parameter of symmetry 
breaking in renormalizable models that utilize the 126, 126 to complete and enforce 
the symmetry breaking down to the SM symmetry |20j. As for Mu it is fine tuned 
to keep a pair of doublets light. The relation between the MSSM couplings at the 
susy breaking scale M s ~ ITeV and the GUT coupling at the scale M x is given by 
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— — = — — — ■ + Sirbiln— f- + AizY^-lnXi - An\i(M x ) (46) 
Oi(Ms) a G {M x ) M s ^ bj 1 v ; 

here 

M° 

X j = l + 8irb j a G (M°)ln-£ (47) 

Ms 

is understood to be evaluated at the values of M x , at G (M x ) determined from the 
one loop calculations. In this equation the contribution of the yukawa couplings has 
not been taken into account and this should also be done in a full investigation |27j . 
Here we will confine ourselves to estimating the corrections using the equations as 
given above , since these were already conjectured[18 to lead to a breakdown of the 
unification scenario. The coefficients 

33 

{&i, &2, & 3 } = (1/16tt 2 ){ — , 1, -3} (48) 

5 

1 ( 199/25 27/5 88/5 \ 

N = 7l6^ 9/5 25 24 (49) 

^ ib7r > \ 11/5 9 14 J 

are the standard one loop and two loop gauge evolution coefficients for the MSSM 
|28j . The term containing Aj represents the leading contribution of the superheavy 
thresholds : 



Xi(fi) = -^-{b lV + b iGB ) + 2{b iV + b iGB )ln— + 2b iS ln— + 2b iF ln— (50) 
21 /i /i /i 

where V,GB,S,F refer to vectors, Goldstone bosons, scalars and fermions respec- 
tively and a sum over heavy mass eigenstates is implicit. The formulae for the thresh- 
old corrections are 



A {th \lnM x ) 
& {th \Log w M x ) 
A {th \sin 2 6 w (M s )) 



5A 1 (M°)+3A 2 (M°)-8A 3 (M° / 



M 
7T 



10&i + 6b 2 - 166 3 
.0217 + .0167(56; + 36' 2 - 8b' 3 )Log 10 j^ 

107ra(M 5 ) 
(56i + 3b 2 - 86 3 ) ^ 3 K Jl 1 X ' 

.00004 - .00024(46; - 9.66' 2 + 5.6b' 3 )Log 10 



(51) 
(52) 



M' 
W x 



(53) 
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Where V i = Wn 2 ^ are the 1-loop beta function coefficients for multiplets with mass 
M' . To evaluate these formulae it is convenient to group the gaugino contributions 
along with the chiral fermions they mix with. The values of the indices Si, 52, S3 
combined as in eqns. fJH^ IHHjl i.e Sw = 4Si — 9.6S2 + 5.6S3; Sx = 5 Si + 3S2 — 8S3 
are given in Table 2 in the appendix. 
The two loop contributions 

A {2 - loop \lnM x ) 
A^- loop \sin 2 e w {M s )) 

Using the values 

a^MxY 1 = 25.6 ; M° x = W 16 - 25 GeV ; M s = ITeV 
a^{M s ) = 57.45 ; oq 1 {M S ) = 30.8 a^(M s ) = 11M (55) 

extrapolated from the global averages of current data, the two loops effects give 

A 2 - l °°v(log w ^f) = -.08 ; A 2 "^( S m 2 ^(M s )) = .0026 (56) 
Ms 

The values of the 1-loop coefficients 6j = lQir 2 bi corresponding to Vector, complex 
scalar and Weyl fermion fields are — llS(i?)/3, S(R)/3, 2S(R)/3 where S(R) is the 
index of the relevant representation. Note in particular that this implies that the 
nonzero hypercharge superheavy vector multiplets which are present in SO (10) models 
will contribute with negative coefficients to the evolution of even the U(1)y coupling. 

We have computed the threshold corrections for a range of values of £ keeping the 
other "insensitive" parameters fixed at randomly chosen representative values 

A = 0.12 ; 7^ = 0.21 ; 7 = 0.23 ; 7 = 0.35 (57) 

The results for different values of these parameters (but with the same £) are 
very similar. We will therefore keep them fixed at these values throughout since here 
we only wish to illustrate the feasibility of precision RG calculations in the SO(10) 
MSGUT. 

For real values of the superpotential parameters the cubic equation (|23|) that 
determines the vevs has one real and two complex (conjugate) solutions. The latter 
give essentially identical corrections. So for real £ we need to present plots for two 
solutions only. These are given as Figs. 1-6. 

From Figs. 1,3 we see that for most real values of £ the threshold effects on 
sin 2 8w{Ms) are less than 10 % of the 1-loop values. There are three exceptional 
values of £ very close to which this limit is breached but even then the change is only 



1 

10&i + 66 2 - 166; 



5&1J + 36 



'2j 



86; 



3j 



InXn 



10na(M s ) b kl 

■7—, ; > £nk(bi — bA—inXi 

(56 x + 36 2 - 863) ^ l3hK 3) k 



(54) 
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Figure 1: Plot of the threshold corrections to Sin 2 6 w vs £ for real £ : real solution 
for x. 

about 25%. For large magnitudes of £ an asymptotic regime of around 10% change 
seems to supervene. Similarly Figs. 2,4 show that the change in Mx is also not 
drastic (though possibly phenomenologically interesting since the gauge contribution 
to the nucleon lifetime goes as M x 4 ) except at certain special points among which 
one recognizes certain known points of enhanced symmetry [1*1 IT3] such as £ = —5, 10 
(SU(5)), £ = 3 (G LR ), £ = -2/3 (flipped S77(5) x 17(1)). It is natural to expect 
that something similar accounts for the other sharp peaks and dips in these plots. 
Moreover their narrowness emphasizes that for generic values of the parameters one 
may expect the threshold corrections to be small for the real £ real x cases. There 
are also regions in which the threshold corrections to Log^Mx are as large as —5 and 
these need special examination with regard to their phenomenological viability and 
consistency with the one scale breaking picture. It is interesting that in this way one 
can scan the parameter space of the MSGUT and obtain a global "tomograph" of the 
variation in its character with the ratio M/m. 

Fig. 3., 4. we give a magnified view of the region |£| < 2. A comparison of the 
graphs shows clearly that the peaks in the threshold corrections coincide by either 
measure, obviously because some particles are becoming very light and enhancing the 
mass ratios that enter the formulae. It will be amusing to use these plots to identify 
and unravel the special regions of the MSGUT parameter space. 
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Figure 2: Plot of the threshold corrected Log^Mx / 'M x vs £ for real £ : real solution 
for x. 

Let us turn next to the complex solutions of the cubic equation for x but still with 
real values of £. We obtain the typical plots Fig. 5,6 The corrections to sin 2 6 w (M s ) 
are very small for small |£| < 2, with a minimum close to £ = 1. From Fig. 6 we see 
that apart from the two peaks near £ = ±5 the corrections to the unification scale 
are quite small for small £ 

When we consider complex values of £ as shown in Fig. 7-12 we see that the 
behaviour is quite regular (like the case of the complex solutions for real £) and once 
again there are large regions of parameter space where the corrections are less than 
10% for Sin 2 9 w while Mx changes by a factor of 10 or less. Thus even this cursory 
scan of the MSGUT parameter space shows that, quite contrary to expectations in 
the literature JH] 5 precision RG analysis of the SO(10) MSGUT is far from being 
futile, since the hierarchy of magnitudes between MSSM one loop gauge coupling 
convergence values (0(a _1 ) effects ) and the one loop threshold and two loop gauge 
coupling corrections (0(1) effects 22 j) is generically maintained at the level of 10% 
or less. Furthermore the RG analysis and parameter scan in terms of the single 
parameter £ can teach us much about the structure of the parameter space since it 
shows a sharp sensitivity to points of enhanced symmetry. We will return to these 
questions at length in the sequel [27] . 
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Figure 3: Magnified Plot of the threshold corrections to Sin 2 9 w vs £ for real £ : real 
solution for x. 



5 d = 5 Operators for B,L violation 
5.1 Higgsino exchange mediated violation 

In SO (10) B — L is a gauge symmetry. Thus the vertices preserve this symmetry and 
the leading effects of the spontaneous violation of B — L by the superheavy S, S vevs 
are just the neutrino mass phenomena. To examine the generation of effective (non- 
renormalizable) operators that violate B+L in the low energy theory via Higgsino and 
gaugino exchange we need the MSSM break-up of the SO(10) invariants 16 ■ 16 ■ 10 
and 16 • 16 • 126. Since we had already presented the PS decomposition of these 
invaraiants in ^l] it is a trivial exercise to use that result to obtain the MSSM wise 
decompositions : 

= V2ti AB [H^T4>B« + H^A^aB ~ + VWA#£fl)] 

= 2\/2h' AB \f, l (tu A d B + Q A L B ) + h(-Q A QB + u A e B - d A 9 B )\ 

- 2V2h! AB hx [d A Q B + e A L B ] + 2^/2h! AB h x [u A Q B + v A L B \ (58) 
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Figure 4: Magnified Plot of the threshold corrected Log\§MxjM\ vs £ for real £ : 
real solution for x. 



™fm = g^ T Cf 7 ii 7* s x£u...* 5 

= 2v / 2(^ a >£x™ - Ej^ft) + 4v / 2Sr A (Cx^ + 

+ <Kt& $ + ^ a \ aX .p) (59) 

whence 



Wj? M = 4:V2f' AB [t 2 (-Q A Q B - u A e B + u A d B ) +t 2 (Q A L B - eu A d B )] 

+ AV2f' AB [-^={h 2 {d A Q B - 3e A L B ) - h 2 {u A Q B - W A L B )} 

+ 2(C 1 d A Q B - C 2 u A Q B ) + 2(E 1 d A L B - D 2 u A L B ) 

+ 2{D 2 e A Q B -E 2 V A Q B )\ 

+ 4f' AB [T l al3 Q Adt Q B + 2i(Ae A e B - G 5 P A P B ) - 2\p2iF x e A v B 

+ (WQ A Q B + 2PQ A L B + V20L A L B ) 

- 2it±(d A v B + u A e B ) + 2iV2(Kd A e B - Jiu A is B )} (60) 
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Figure 5: Plot of the threshold corrections to Sin 2 9 w vs £ for real £ : complex solution 
for x. 



We have suppressed G321 indices and used a sub multiplet naming convention 
specified in Section 2. and and Table I in the Appendix. 

In order that the exchange of a Higgsino that couples to matter with a given 
B + L lead to a B + L violating d = 5 operator in the effective theory at sub GUT 
energies it is necessary that it have a nonzero contraction with a conjugate (MSSM) 
representation Higgsino that couples to a matter chiral bilinear with & B + L different 
from the conjugate of the first B + L value. Inspecting the above superpotentials 
one finds that only {f(i),?(2)} and {£(1), £(2), *(4)} satisfy this requirement. Terms 
containing the right handed neutrinos V\ must be further processed to integrate out 
the heavy field Pa m favour of the the light neutrinos v& . This will introduce an extra 
factor of rn v Dira J M v Ma j orana and effectively lead to amplitudes suppressed like those of 
d = 6 operators. Thus on integrating out the heavy triplet Higgs supermultiplets one 
obtains the effective d = 4 Superpotential for Baryon Number violating processes : 

W tff^ = ° = l abcd(^QaQbQcL d ) + RABCD^e A u B u c d D ) (61) 
where the coefficients are 
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Figure 6: Plot of the threshold corrected Log^Mx / 'M x vs £ for real £ : complex 
solution for x. 

Labcd = S^fiAshcD + S\ 2 hABfcD + S 2 l fAshco + <S 2 2 fAsfcD (62) 

and 

Rabcd = S x Kab^cd ~ S\hAsfcD — S 2 X f ab^cd + S 2 'J ? abIcd 

- iV^S^fABhcD + iV^S^fABfcD (63) 

here S = T -1 and T is the mass matrix for [3, 1, ±2/3]-sector triplets : W = 
iTt + . . . , while 

h AB = 2V2h' AB f ab = 4V2f AB (64) 

This expression and the "Clebsches" contained in it , as well as the new baryon 
decay channel mediated by the triplets (t( 4 )) contained in £ 12 6(10, 1,3) ( the same 
PS multiplet that contains the Higgs field responsible for the right handed neutrino 
Majorana mass), were given in jllj . Previous work ^2] on S126 mediated decay 
focussed on the multiplets v?' ,tP> and found that there was no contribution of t^\t^ 
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Figure 7: Plot of the threshold corrections to Sin 2 9 w vs Re(£) for complex £ : 
Im£ — 1.2 , first solution for x. 



in their models. This new channel nominally strengthens the emergent link between 
neutrino mass and baryon decay. Note however that couples only to the RR 
combinations dv + ue and as such its exchange will contribute only to the RRRR 
channel which , at least in SO(10), seems [12] generically suppressed except at very 
large tan/3. However the mixing in the triplet mass matrix could also strengthen the 
effects of this channel. 

5.2 Novel d = 5, A(B J r L) ^ Operators via superheavy gaug- 
ino exchange ? 

A novel situation apparently arises in this GUT due to exchange of superheavy gaug- 
ino Dirac multiplets that couple to matter both via the gauge yukawa couplings of 
their gaugino part and the superpotential couplings of their (126) chiral components 
to the matter sector. Such gauginos are not present in the case of SU(5) As is evi- 
dent from eqn.flUl) the 126 submultiplet fields E 2 [3, 2, -1/3], I, -2], J x [3, 1, 4/3] 
(which mix with the superheavy SO( 10) /SU(5) coset gauginos : see Section 2. and 
3.) couple only to terms containing at least one superheavy neutrino field va- Thus , 
to leading order in My 1 , the exchange of such gaugino dirac multiplets will not lead 
to any d = 5 operator with 4 light external fields. However a puzzle remains. 
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Figure 8: Plot of the threshold corrected Log\§MxjM\ vs Re(£) for complex £ : 
7m£ = 1.2 , first solution for x. 



The superheavy neutrinos mix with the usual light neutrinos via Dirac masses. 
So in the effective theory one trades them for the light neutrinos by using their 
equations of motion to leading order in their (inverse) masses ( effectively v = 

-2( m Dirac/ M Majorana^ + y partg ^ J x of the Gauge DhaC E,F,J 

multiplets therefore couple to light neutrinos and another light matter field with a 
small coupling ~ 0{m Dtrac / M Ma j orana ) . Exchange of the gaugino dirac fermion be- 
tween a gauge yukawa vertex and a 126 -16-16 vertex can lead to effective operators 
involving 4 light matter fields of which at least one is a light neutrino and one is 
anti-chiral. This appears to violate the usual argument that in the effective MSSM 
arising from a Susy GUT , supersymmetric D terms involving 4 light (mixed chiral and 
anti-chiral) fields must be d > 6 or equivalently that the d = 5, B, L violating terms 
are either of form [QQQL]p or [euud}F- Exchange of S'O(IO)/ SU{h) coset gauginos 
peculiar to 5*0(10) however appears to lead to (admittedly suppressed) d = 5 chiral- 
anti-chiral operators with 4 light fields. These operators arise once the Electroweak 
scale vev that gives rise to neutrino Dirac masses is turned on. This vev is smaller 
than M s and arises after soft susy breaking terms are included. In this theory B — L 
is spontaneously broken giving rise to the Majorana mass for conjugate neutrinos 
(which was used to eliminate them in favour of the SM neutrinos). Thus perhaps the 
contradiction is not as violent as it seems at first. We emphasize that there are no 
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Figure 9: Plot of the threshold corrections to Sin 2 9 w vs Re(£) for complex £ : 
Im£ = 1.2 , Second solution for x. 

analogous processes in SU(5) Susy GUTs since there the 12 coset gauginos acquire 
partners from the purely AM type 24 plets which do not couple to the matter sector. 

The couplings of the gauginos of SO (10) to the matter fields are easily computed 
by adapting the PS reduction of the SO(10) covariant derivative for the spinor 16 

C g .y = 2^[^(^)^a+^A A ((-^), K )^% 

+ ^[tA^^ Q +CA;/ ft ]+ff.c (65) 
The terms carrying B, L are are 

£a(b+l)^o — \/2g[L* J 4 Q + Q*J 4 L] — V2g[d J 4 e + u J 4 z/ + e* J 4 d + v* J±u\ (66) 
+ Q?/V2) [- d*X 3 a L - uE (b) L + TX 3 Q + d*E 5 Q + ed E 5 Q + eu X 3 Q] 
+ (g/V2)[(L*(X 3 d- E 5 u) - Q*(X 3 e - E 5 u) + eQ*(E 5 d- X 3 u)} + .... 

(67) 
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Figure 10: Plot of the threshold corrected Log^Mx / 'M^ vs Re(£) for complex £ : 
Im£ = 1.2 , Second solution for x. 



There are no X[3, l,±5/3] sector submultiplets in the 126. Thus we can focus 
on just the E[3, 2, ±1/3] and the J[3, l,±4/3] sectors here. As discussed in Sec- 
tion 2., the superheavy gauginos J^,E 5 mix with 126 derived fermions Ji [3, 1,4/3] 
and £2 [3, 2, — 1/3]. Examining eqn lffiUJ) we see that Ji,E 2 couple only to operators 
involving at least one superheavy v field [E\ G 126 does not mix with E-gauginos): 

W§ M = -8V2 f' AB [E 2 Q B + %J x u B \v A = [E 2 Qa + iJiu A }{f'M^m vD ) AB v B + ... (68) 

Since J\ couples to a B = — 1/3, L = 1 operator while J4 couples only to B — 
1/3, L = —1, J exchange does not lead to B + L violation. The E sector gaugino i.e 
E 5 couples as 

-^=[e~dE 5 Q + Q*E 5 u-L*E 5 u] (69) 
v2 

Only the first terms in ()68I69|) are relevant and thus we find the following effective 
lagrangian due to superheavy gaugino exchange : 

L A (b + l)=2 = ±gU\M v )- l m^) AB £-\ 2 [ed c Q c {Q A v B + Q A v B )\ (70) 
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Figure 11: Plot of the threshold corrections to Sin 2 9 w vs Re(£) for complex £ : 
7m£ = 1.2 , Third solution for x. 

where is essentially the mass of the exchanged gaugino times mixing factors 
written compactly interms of the inverse of the relevant fermion mass matrix (in the 
E[3,2,±l\ sector). By dressing this with MSSM gauginos we obtain AB = AL = 1 
violating 4 fermi vertices responsible for processes like 

uidL — >d L + T7Z (71) 

This is a vertex quite distinct from the Higgsino mediated vertices since it involves 
exchange of massive gauginos between a chiral and an anti-chiral vertex. It requires 
non zero external momenta for the fermions and vanishes in the limit of zero external 
momenta. Thus the coefficient of the corresponding 4-fermi operators for B violation 
in the effective lagrangian is ~ M vL } m Nud / 'M\M\ where M s is the Susy breaking 
scale. This magnitude seems hopelessly suppressed (relative even to gauge boson 
exchange) to be observable. Nevertheless the contrast of its structure with that of the 
standard QQQL and uude operators perhaps warrants a more thorough investigation 
of the conditions for the possibility of its appearance in the effective theory. 
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Figure 12: Plot of the threshold corrected LogioMx/Mx vs Re(£) for complex £ : 
im£ = 1.2 , Third solution for x. 



6 Fermion Mass Formulae 

A vital issue for any SO(10) GUT is the type of predictions it makes for the relations 
among the parameters of the (Type I and Type II) seesaw mechanisms [16J by which 
Neutrino masses and mixings arise. From the coupling of neutrinos to the 126 we 
find that the Majorana mass matrix of the superheavy neutrinos i>a is ( eqn. (j60|0 

M P AB = -4iV2f AB < >= iV2f AB W (72) 

Similarly the Majorana mass matrix for the left neutrinos v A is ( eqn. (|60Jl ). 

M» AB = 4V2f' AB < O u >= 8if AB <0 > (73) 

where < 0_ > is the small vev of the SU{2) L triplet in the (10, 3, l)^ induced by 
a tadpole that arises as a consequence of SU(2)l breaking(see below). 

In addition to this there is the Dirac mass which mixes the left and right neutrinos 

m% = 2V2ti AB < > +4iV6f AB < hf > (74) 
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We must make the fine tuning DetTi = necessary to keep a pair of Higgs dou- 
blets Hm,H(ij (which is to develop the EW scale vev ) light. Then these doublets 
will be the left and right null eigenstates of the mass matrix 7i. If the bi-unitary 
transformations responsible for diagonalizing and TiH) are U, U i.e 

Diag(m { i\m%\ ....) = U T HU (75) 

then writing 

h {i) = U l3 H ij) ; fc« = UijH® (76) 

where H®,H(f> are the mass eigenstate doublets, the contributions of any cou- 
pling in which h^' enter can be accounted for in the effective MSSM below the 
heavy thresholds of the GUT just by replacing h^' — > otiH^ ,h^' — *> oiiH^ 1 ' where 
oci = Un,C(i = Un and we have numbered the massless doublet pair "1". These com- 
ponents are easily obtained from the normalized null eigenvectors V, V of 1-0 Ti and 
HTO to be Un = Vi, Ua = V*. Thus the neutrino Dirac mass matrix becomes 

Mf B = {2V2h! ABai + iiV6f' AB a 2 )v u (77) 

To obtain the final formula for the neutrino masses and mixings we must eliminate 
the v fields which are superheavy and evaluate the tadpole that gives rise to the Type 
II seesaw. The first step is standard. As for the O(10, 3, l)i 2 6, O(10, 3, l)^ vevs, 
inspection of the mass spectrum (Table I in Appendix) and eqns. ()36ll38|) yields the 
relevant terms in the superpotential as 



pm = M d-d-^H a ^ 4a O^-^=H aa ^O^ 



since £ 4 4ai = ^h^^ 2 = ^/if 2) ,$ 44iQ = -V2h£\$£ = y/2h£\ one gets for 
the relevant terms : 



ol^v\ — + V2(i'yai + i2v / 3^a 2 )a 4 f^ (79) 

Thus the vev we need i.e < CL > is immediately determined to leading order in 
Mw/Mu by by the equation for 0+ as 

v 2 

< 0_ >= V2(»7ai + 2iV3ria 2 )a 4 {-^) (80) 

Mo 

and Mo can be read off from Table I to be Mo = 2(M + r](3a — p)). 
The quark and charged lepton mass matrices are 
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M 



M d 



M' 



u 




(81) 



(82) 



(83) 
(84) 



These formulae are now in a form ready to use for fitting the fermion mass and mixing 
data after lifting it via the RG equations of the MSSM to the GUT scale. 

7 Discussion and Outlook 

In this paper we have calculated the complete superheavy spectrum of the Mini- 
mal supersymmetric GUT along with the gauge and chiral couplings of all MSSM 
multiplets in a readily accessible form. Partial calculations of these spectra and cou- 
plings [31 El El E] have been published earlier but our method is different from 
the computer based method of [3J El E] and is more complete, especially regarding 
couplings. Being analytic and explicit it also allows us to trace and resolve discrep- 
ancies arising within the computer based approach. We used the calculated spectra 
to perform a preliminary scan of the parameter space of the MSGUT as regards the 
magnitude of the threshold corrections to two crucial phenomenological parameters 
of the MSGUT : the Weinberg angle at low energy and the mass of the X lepto-quark 
gauge supermultiplet. We obtained a result that is in sharp contrast with expecta- 
tions in the literature^H] that precision RG calculations in SO(10) are futile. On the 
contrary we find that the 1-loop GUT threshold and gauge two loop contributions 
are modest but significant. Thus, on the one hand, the basic GUT picture suggested 
by the convergence of gauge couplings in the MSSM is in fact not destroyed by the 
contributions of the large number of superheavy fields. On the other hand extant 
precision calculations that ignore threshold effects in SO (10) GUTs seem to be of 
dubious validity. In particular the proper RG analysis of the MSGUT taking into 
account EWRSB, all fermion masses and GUT threshold effects still remains to be 
done. This calculation is now being performed [27\. In view of the other phenomeno- 
logical successes of renormalizable Susy SO(10) [SHEIE! and the unforseen correlations 
between disparate phenomena like neutrino oscillations and nucleon decay that have 
emerged |12[ ITT] the mildness and calculability of threshold effects in the MSGUT is 
a most welcome and promising development. Our preliminary scans of the MSGUT 
parameter space (whose very feasibility - based on there being just one "sensitive" 
control parameter (£) - is a matter of some astonishment) show that the threshold ef- 
fects can potentially narrow down the allowed regions of the MSGUT parameter space 



35 



and indicate correlations between the GUT scale and the B — L violating scale which 
can be of crucial significance when cross checking the particle physics phenomenology 
against cosmology. We have also argued that above the perturbative unification scale 
realistic renormalizable SO(10) GUTs are necessarily strongly coupled [HUH]- We have 
recently reported the results of 2-loop calculations of MSGUT RG equations above 
the SO(10) restoration scale which we used to show that the strong growth of the 
SO (10) coupling above Mx cannot be evaded by taking shelter in a weakly coupled 
fixed point [201 • On the other hand our work [8, 9. has shown that a scenario of a cal- 
culable dynamical symmetry breaking of the GUT symmetry which utilizes the nearly 
exact supersymmetry at the GUT scale offers rich possibilities for the significance of 
the new length scale associated with the condensation of SO(10)/Gs2i coset gaugi- 
nos implied by both holomorphic analysis and by the Konishi anomaly. The present 
calculations show the way for crossing the threshold and entering into the SO (10) 
regime in a controlled way. The emerging coherence of the low energy phenomenol- 
ogy, B and L violation, perturbative GUT structures (such as the natural R-parity 
preservation jSU 120] , successful seesaw scenarios, leptogenisis etc) and exciting hints 
of deeper mysteries, perhaps unveilable[Hl E], carry to our hopeful nostrils the spoor 
of a grail perhaps within reach. 



Note Added 



After this paper was posted on the arXive as hep-ph/0405074 the authors of 
|hep-ph/0 405300 claimed that the mass spectra listed in Appendix A were not inter- 
nally consistent with the requirements of SU(5) or SU(5) x U(l) symmetry (at the 
special vevs where p = a = ±u). However this is incorrect. The mass spectra we 
derived via a PS decomposition of SO(10) organize straightforwardly and termwise 
into appropriate SU(5) invariants for SU(5) invariant vevs as given in the Appendix 
B (added as above to hep-ph/0205074vl). This term- wise reorganization of several 
hundred G123 invariant mass terms into SU(5) invariant mass terms is a more strin- 
gent consistency test than the tests of hep-ph/0405300 which are based on traces and 
determinants and valid only for their conventions. The phase conventions and field 
normalizations of hep-ph/0405300 are quite different from our work. Thus the blind 
application of their trace and determinant consistency tests to our results cannot but 
fail. We maintain unit field normalizations throughout by using only unitary field re- 
definitions of fields with canonical kinetic terms . Finally our results for chiral spectra 
also coincide, up to minor convention related adjustments, with those obtained in a 
parallel computation reported in [13] . 

Finally we stress that our method [TT] yields all coupling coefficients between both 
spinor and tensor irreps and not just the tensor irrep ones relevant for masses and 
symmetry breaking which were obtained using the method of JU]. Moreover we 
note that the most complex of the mass matrices given here namely those of the 
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Higgs doublets and triplets relevant for proton decay were already derived by us in 
|hep-ph/0204097> 2(2003) P] 

Further Note Added 

After version 2 of this paper (including the SU(5) reorganization given in Appendix 
B above) was accepted for publication the authors of hep-ph/0405030D issued yet 
another preprint(hep-ph/0412348 vl ) this time claiming that although our results 
pass the SU(5) reorganization 'test ' for a = o = they failed to do so for a = a ^ 
and that the counting of Goldstone modes and distinct mass eigenvalues was, in their 
opinion incorrect. Further they claimed that our results were inconsistent since they 
failed to pass certain 'trace and hermiticity tests ' that they had applied successfully 
to their own results. All these claims are incorrect. Our results in fact pass all three 
tests. We have issued a preprint showing this explcitly [33]. Here we only remark that 
once the super-Higgs effect for 5*0(10) — ► MSSM has been verified it is scarcely 
feasible that the Goldstone- Higgs counting could fail for SO (10) — > 577(5) since the 
latter is a special case of the same spectra ! However the reader can easily check that 
the SU(5) singlet and lOplet mass matrices have zero determinant confirming that the 
required Goldstone supermultiplets 1 + 10 + 10 are present. The demonstration that 
the trace constraints and 'hermiticity' tests of hep-ph/0205300 are also satisfied is also 
straightforward once proper account is taken of the difference in the phase conventions 
of the two calculations. Details may be found in [33J. Finally as this paper goes to 
press the authors of hep-ph/0412348vl have reissued the preprint hep-ph/0412348v2 
in which all the claims of the inconsistency of our results are totally retracted. 
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Appendix A : Tables of masses and mixings 

In this appendix we collect our results for the chiral fermion/gaugino states, masses 
and mixing matrices for the reader's convenience. Apart from the discussion of gauge 
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multiplet masses our results have been obtained in parallel with and are compatible 
with those of which, however , are computed with a different normalization 
for the E, E fields resulting in a difference between the mass and yukawa coupling 
parameters M, r\ of these multiplets in the two starting actions. Moreover certain 
minor phase differences also exist between the definitions of representative states 
used by them and our definitions for the same states (which follow directly from our 
consistent definitions of PS tensors from SO(10) submultiplets). Mixing matrix rows 
are labelled by barred irreps and columns by unbarred. Unmixed cases(i)) are given 
as Table I. 

ii) Chiral Mixed states 



a) [8, 1,0] (i^ifc 



'p. > ( Pp,(R0)J 



TZ = 2( ( m T-^ ~^ XU A (85) 
V -V2\u m + Xip-a) J K ' 



m R± = \n ± \ = \2m[l + (| - a) ± y (|) 2 + 2u>*]\ (86) 
The corresponding eigenvectors can be found by diagonalizing the matrix VSR) . 
b) [1, 2, -1] (h lt h 2 , h 3 , h A ) ©[1,2, 1] [hx, h 2 , h 3 , h A ) 

— fxra V (15)Q V( 15 ) Q feffifff T {15) V^ 15 ) $2L\ 



n 



—Mh +7v / 3(w — a) — 7V / 3(cj + a) — 70" ' 

-f\/3(w + a) -(2M + 4r](a + u)) 

7v%(cu-a) -(2M + 4r](a - u)) -2rfo^fl 

\ -0-7 -2t)o*J% -2m + QX(uj-a) J 

The above matrix is to be diagonalized after imposing the fine tuning condition 
DetTi = to keep one pair of doublets light. 



C ) 1) fjK^b ^2, ^3? ^4, £5) © [3, 1, — jjK^Ij ^3) ^4, £5) 

= (-^ A4 5 ^fo)' S fo)> S fiO' ( f ) i(R+)) © (^M> S (a)/i4, E^ 4 ( a ) , E^ 4 ( fl o) , <f>ftR-)) 



( 




7(a + p) 7(p-o) 


2\/2iu*/ 


zo"7 






7(P - a ) 


2M 












7(p + a) 


2M 


Ay/2iur) 


2irfa 






-2V2za;7 


-AV2iur] 


2M + 2r/p + 2r/a 


-2V2rja 




V 


zo"7 


2zr/a 


2^(7 


-2m - 2X(a + p 


-4u) ) 



iii) Mixed gauge chiral . 
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Fie\d[SU(3),SU(2),Y] 



PS Fields 



Mass 



A[l,l,4],A[l,l,-4] 
d[8, 2,1], (74(8, 2,-1] 

C 2 [8,2,l],(7 2 [(8,2,-1] 
D 1 [3,2,l\,D 1 [(3,2,-l] 

D 2 [3, 2, J], L> 2 [(3,2,-|] 

^[3,2, |],^[(3,2, — |] 
K[3,l,-f],J?[(3,l,§] 

L[6,1,|],Z[(6,1,-|] 



M[6,1,|],M[(6,1,-|] 
7V[6,1,-|],7V[(6,1,|] 

0[l,3,-2],0[(l,3,+2] 
P[3,3,-|],P[3,3,|] 

^[6,3,|],W[(6,3,-|] 

7[3,l,f],7[(3,l,-f] 

5[1,3,0] 

Q[8,3,0] 

C/[3,3,|],C7[3,3,-|] 



^[l,2,-3],t/[l,2,3] 
S[6, 2, |], ^[(6,2,-1] 

r[6,2,-|],y[(B,2,|] 

Z[8,l,2],Z[8,l,-2] 



f%R+) ^44(fl-) 

v^2 ' V2 
y A y A . 

y 1 . y A 

y ■ y v 



y 4 y 



y ■ y v 

/v=r'(-R0) \ 

l^P 5 ^(RO))^ 

y'__ — s mm 

/vt'(-R+) v'^" 5 ^ 
\^fiV(R+)i ^(R-)JP<i 

/vt'(-R — ) v-i'/iv ^ 
K^ixv ■> ^(R+))p.<v 



/IV 

_2.4 

^44(L) S (L) 



^pA(L), ^(L) 



p,D(L)i^(L) 
$ P(R+)' ^4(i?-) 
<?(L) 

0p(L) 5 4 (L) 



^44g2 ^gj 

V2 ' v / 2 ; __ 
<t> V n {R+)<PrtR-) 



2(M + r](p + 3a + 6u)) 
2(-M + r](a + uj)) 

2(-M + r)(a-u)) 
2(M + r)(a + uj)) 

2(M + 7](a + 3a;)) 

-2(M + 77(0-0;)) 
2(M + ?7(a + p + 2cj)) 

2(M + ??(jo-a)) 



2(M + ?7(p-a + 2cj)) 
2(M + r?(p - a - 2a;)) 

2(M + r/(3a-p)) 
2(M + r/(a-p)) 

2(M -r](a + p)) 
-2(m + \(p + a + 4lu)) 

2{m + \{2a - p)) 
2(m - X(a+p)) 
—2{m — \{p — a)) 

2(m + 3A(a + uo)) 
-2{m + A(u; - a)) 
2(m - A(a + a;)) 
2(m + A(p - a)) 



Table 1: i) Masses of the unmixed states in terms of the superheavy vevs . The 
SU(2) L contraction order is always F a F a . The primed fields defined for SU(3) C 
sextets maintain unit norm. The absolute value of the expressions in the column 
"Mass" is understood. 
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FieldfS'f/f3 N ) SU(2) Y] 


'-'2) u l J 




^ X 


AFl 1 41 


Jn 12/51 

l U > u, ±z/ o j- 


9 6 


12 


ftlfi 2 5/31 


15 3 51 


19.2 


—6 


(7 8 2 1 


16 4 12/51 


4 8 


—24 


D\?> 2 7/31 


11 3/2 49/101 


10 8 


21 


M3 2 1 /3l 


11 3/2 1 /10l 

\L, O/ Z, 1/ 1UJ- 


—8 4 


■J 


F 1 1 21 


10 3/51 


2.4 




G\l 1 

1, UJ 


10 01 

|U, U, U J 


n 


n 

W 


Ml 2 11 

/t[±, Z, 1J 


10 1 /2 3/101 


—3 6 


3 

*J 


7[3 1 10/31 


11/2 51 


22.8 


21 


7T3 1 4/31 


11 /2 4/51 


6 


n 


M3 1 8/31 

-fY [O, 1, O/ OJ 


11/20 16/51 
\1/ z, u, iu/ oj- 


15 6 


12 


fifi 1 2/31 


15/2 2/51 
\<J/Z, u, z/oj 


1 5 6 

X tj . VJ 


— 18 

X(J 


M[6 1 8/31 


15/2 32/51 


39 6 


12 


A/Tfi 1 4/31 


15/2 8/51 
|0/ z, u, o/ 


20 4 

_ VJ . t: 


— 12 


OTl, 3, 21 


10,2,9/51 

L ' ? / J 


-12 


15 


P[3,3,2/3] 


{3/2,6,3/5} 


-46.8 


9 


Q[8,3, 0] 


{9,16,0} 


-103.2 


-24 


#[8,1,0] 


{3,0,0} 


16.8 


-24 


5(1,3,0] 


{0,2,0} 


-19.2 


6 


t[3, 1,2/3] 


{1/2,0,1/5} 


3.6 


-3 


C/[3,3,4/3] 


{3/2,6,12/5} 


-39.6 


18 


\/[l,2,3] 


{0,1/2,27/10} 


6 


15 


W[6,3,2/3] 


{15/2,12,6/5} 


-68.4 


-18 


X[3,2,5/3] 


{1,3/2,5/2} 


1.2 


9 


F[6,2,l/3] 


{5,3,1/5} 





-30 


Z[8,l,2] 


{3,0,24/5} 
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Table 2: Index values for the 26 different chiral multiplet types (used in the threshold 
corrections). Except for Q,R,S all other reps come in complex pairs. Sw = 4 Si — 
9. 6S2 + 5. 6S3, Sx = 5 Si + 3S*2 — 8S3 are the combinations that enter the threshold 
corrections to Sin 2 9 w and to Log w M x 
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a) [1, 1, OKd, G 2 , G 3 , G 4 , G 5 , G 6 ) ee (0, 0( 15 \ 0^, %, ^f±) , VS^WS*") ) 



£ = 2 



/ m 





v^6Ac<j 


V2 




\ 





m + 2Aa 


2y / 2Au; 


—ir]^/3a 


—irja^ 







2-\/2Au; 


m + A(p + 2a) 





irja 

7f 




— ir)\/3a 





M + r](p + 3a- Quj) 


V5gcr* 
2 




—ir\a^\ 


ir)y/3a 


M + ??(j9 + 3a - Quj) 





2 


V o 








VEga* 
2 


\flga* 
2 


o / 



b) [3, 2, -\\(E 2 , E 3 , E 4 , E 5 ) © [3, 2, §](£ 2 , £ 3 , £ 4 , E 5 ) 

= ( S 4ai ? 0[l 4 ) a2 ' ^S' i4 ' A 3 © ( S Ja2' ^4ai' ^jSai' "W) 



£ = 



^ -2(M + r?(a-3cj)) 
2iV2r/a 
—2irfa 



-2\j2irio 
-2{m + A(a - to)) 
-2y/2Xv 
2g(a* - uj*) 



2ir]<T 
-2^/2\uj 
—2{m — \uS) 
gV2(u* - p*) 



c) [1, 1, -2KA, F 2 , F 3 ) © [1, 1, 2] (Fx, F 2 , F 3 ) 

= (S44(fl0), © ( S (J?0)> • 



2<?(a*-co*) 




(87) 



/ 



/ 2(M + + 3a)) -2iy/Zrp -gV2a* \ 
JF = 2iv%j7a 2(m + A(p + 2a)) v 7 ^^*) (88) 

V ~gV2a* -V24iguj* / 

d) [3, 1, -|](Ji, J 2 , J 3 , J 4 ) © [3, 1, |](Ji, J 2 , J 3 , -A) 

= 4 MM) , A/) © (£ W) , 0% 0^ o) , AJ) 



2(M + r](a + p-2cj)) -2rp 

2r](T —2(m + Aa) 

-2y/2rja -2^2\uj 
\ -igV2a* 2^2iga* 



2V2rja -igV2a* 

-2^/2\uj -2igV2a* 

-2{m + A(a + p)) -Aiguj* 

Aiguj* 



e) [3, 2, §] (X u X 2 , X 3 ) © [3, 2, - §] (X U X 2 , X 3 ) 



>' ^/J4a2' A p4a2J 



?9) 



2(m + A(a + cu)) -2y/2\u -2g(a* + uo*) 
X = I -2v^Acj 2(m + Acu) >/2ff(w*+P*) 

-2#(a* + cj*) v^^+P*) 



(90) 
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Appendix B : SU(5) x U(l) Reassembly Crosscheck 

Given the complexity of the spectra and couplings derived here it would be useful 
to have a method of cross checking the internal consistency of our results. A stringent 
check is provided by verifying that at special values of the vevs i.e 



where the unbroken symmetry includes SU(5) the MSSM labelled mass spectra and 
couplings given in Appendix A do indeed reassemble into SU(5) invariant form. For 
the mass spectra this is fairly straightforward to check and is reported explicitly below. 
A similar calculation [22] for the super potential couplings is much more tedious but 
furnishes an SO (10) - SU(5) x U(l) analog of the "SO(10)-PS Clebsches " reported 



The decomposition of the chiral multiplets of the MSGUT into SU(5) x £7(1) 
multiplets and of those into MSSM multiplets (named as per the alphabetic conven- 
tion of Appendix A) is given below. The only complication is that certain MSSM 
multiplet types occur in several copies and (orthogonal) mixtures of these are present 
in the different SU(5) mutiplets. Thus, for instance, the 210 contains a 24 and a 
75 of SU(5) both of which contain mixtures of the G123 multiplets Ri(8, 1,0) and 
#2(8, 1,0). These mixtures must be orthogonal and must be precisely the eigenstates 
of the mass matrices in this G123 sector which have the same masses as the rest of the 
G123 submultiplet sets within the 24-plet and 75-plets as two wholes. The fact that 
this follows in every case from our results appears to confirm their reliability. The 
decompositions we need are : 



p = a = ±cu 



(91) 



here . 



H 




5-i 



10-3 



15 3 



45x 



50_i 



-)+C 2 (8,2,-l) 



S 



= 126 = 1 5 (G 5 ) + 5i + 10 3 + 15_ 3 + 45_i + 50! 
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5! = /i 2 (l,2,l)+* 2i4 (3,l,--) 

T0 3 = A(l, 1,-2) + Ji(3,l,^) + ^ 2 (3,2,-^) 

15-3 = 0(1, 3, 2) + ^(3,2,^ + ^(6,1,-^) 

45_! = ^(i,2,-l)+t 2 (3,lJ) + P(3,3j) + X(3,l,-^)+ J D 1 (3,2,^) 

+ L(6, 1,^+^(8, 2,-1) 

50x = A(l, 1, -4) + t 2i4 (3, 1, ~) + L» 2 (3, 2, -I) + W(6, 3, - H) + M(6, 1, jj) + C 2 (8, 2, 1) 
$ = 210 = 1 + 5_ 4 + 5 4 + 10 2 + T0_ 2 + 24 + 40 2 + 40_ 2 + 75 

lo — Oi j2j 3 

5_ 4 = / i4 (l,2,l) + t 5 (3,l,-^) 

5 4 = h4(l,2,-l) +£5(3,1,3) 

10 2 = F 2 (l, 1, 2) + J 2)3 (3, 1,-^+^3,4(3, 2, 1) 

T0_ 2 = F 2 (l, 1,-2) + ^3(3,1,^+^3,4(3, 2,-^) 

24 = (1, 1, 0)Gi )2)3 + 5(1, 3, 0) + X 1>2 {3, 2, -|) + X 1>2 (3, 2, jj) + i2 1)2 (8, 1, 0) 

40 2 = V(l, 2, -3) + £3,4(3, 2, 1) + J 2 , 3 (3, 1, ~) + C7(3, 3, ~) + Z(8, 1, 2) + F(6, 2, 1) 

40_ 2 = t/(l, 2, 3) + £3,4(3, 2, -i) + J 2 , 3 (3, 1, ^) + 17(3, 3, ^) + Z(8, 1, -2) + F (6, 2, -1) 

10 - - 10 5 - - 5 

75 = (l,l,0)G li2i 3 + /(3,l, y ) + /(3,l,- y )+X li2 (3,2,--)+X li2 (3,2,-) 

+ 5(6,2,^) + J B(6,2,-^) + J R 1 , 2 (8,l,0) + g(8,3,0) (92) 

If we insert a = — uj = p in the mass matrices of Appendix A we find that, after 
diagonalizing the mass matrices of the submultiplets that mix, the resultant spectra 
group precisely as indicated by the decompositions above with all the subreps of a 
given SU(5) irrep obtaining the same mass. One obtains the SU (5) invariant mass 
terms : 

2(M + Wrip)!^ + 2(M + Ar]p)5 s 5^ + 2(M - 2?7p)50 E 50s 
+ 2(M + 4?7p)10 s l% + 2(M + 2?7p)l5 E 15s + 2M45 S 4% 
+ M H 5 H 5 H + (m + 6Ap)(l$) 2 + 2(m + 6Ap)5*5* + 2(m + 3 Xp) 10 $TU$ 
+ (m + Xp)(24$) 2 + 2m40< E> 40$ + (m - 2Ap)(75$) 2 
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+ 2r/v / 3(cf(5 s 5$ + 10 s 10$) + ^(5^5$ + lOglO*)) (93) 
+ 2v / 3p(75 s 5h + j5 H 5^) + 2zry\/5(crl^l$ - + ^ya5^5 H + W^hS® 

Where every 577(5) invariant has been normalized so that the individual G123 
sub-rep masses can be read off directly from the coefficient of the invariant for com- 
plex SU(5) representations which pair into Dirac supermultiplets and is 2 times the 
coefficient for the real representations which remain unpaired Majorana/Chiral su- 
permultiplets. For a = —uj = p,cr = a = the 20 Goldstone supermultiplets 
G,J,J,F,F,E,E) of the coset SO (10)/ 577(5) x U(l) remain heavy (see the gauge- 
chiral super-Higgs mass formulae in Section 2) as they should since they are eaten 
in the spontaneous breaking SO(10) — ► SU(5) x U(l) while the 12 fields in the 
{X3,X3} multiplets lose their mass terms with {Xi^, X\^} since they form part of 
the unbroken SU(5) gauge supermultiplet. When a = u = p i.e for flipped SU(5), the 
roles of the {X 3 ,X 3 } and {E 5 ,E 5 } gauge multiplets are interchanged, with the E's 
remaining massless and the X's becoming heavy, so that one obtains the SU(5) in- 
variant groupings corresponding to the "flipped" SU(5) xU(l) C 5*0(10) embedding. 
Note that this successful SU(5) reassembly is a much more fine-grained consistency 
test than any overall trace or determinant test. 
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